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A large number of recent works point to the emergence of intriguing analogs of fractional quantum Hall 
states in lattice models due to effective interactions in nearly flat bands with Chern number C = 1. Here, 
we provide an intuitive and efficient construction of almost dispersionless bands with higher Chern numbers. 
Inspired by the physics of quantum Hall multilayers and pyrochlore-based transition-metal oxides, we study a 
tight-binding model describing spin-orbit coupled electrons in N parallel kagome layers connected by apical 
sites forming N — 1 intermediate triangular layers (as in the pyrochlore lattice). For each N, we find finite 
regions in parameter space giving a virtually flat band with C — N . We analytically express the states within 
these topological bands in terms of single-layer states and thereby explicitly demonstrate that the C = wave 
functions have an appealing structure in which layer index and translations in reciprocal space are intricately 
coupled. This provides a promising arena for new collective states of matter. 

PACS numbers: 73.43.Cd, 71.10.Fd, 73.21.Ac 



Introduction. Topological insulators, predicted 

theoretically^ a few years ago and subsequently experi- 
mentally observed,^ have attracted enormous amounts of 
interest. As these systems can usually be understood in a 
single-particle picture, an intriguing question is whether 
interactions can lead to qualitative new phenomena. Evidence 
that this is indeed the case has been accumulating during the 
past year. In a key step it was shown that e.g., an appropriate 
combination of ferromagnetism and spin-orbit coupling can 
lead to nearly flat bands with unit Chern number, C = 1, 
in itinerant lattice sy stems. Given the flat bands, these 
systems are likely to host interesting strongly correlated states 
and, at least in theory, this opens up a number of intriguing 
perspectives including high-temperature fractional quantum 
Hall states^ . Indeed, numerical exact diagonalization studies 
convincingly show the existence of such states. While the 
list of flat band models with C = 1 is still growing and 
a better understanding of the relation between these bands 
and Landau levels^' is developing, bands with higher 
Chern number could host qualitatively new phases of matter^^ 
as they have no direct analogue in the continuum.^^ In this 
context, two very recent papers made interesting progress 
by introducing models harboring relatively flat bands with 
C = 2.24'25 

At the same time there is a tremendous experimental de- 
velopment and theoretical interest in the physics of transition 
metal oxides, especially the iridates, where the effect of spin- 
orbit coupling is profound.^^'^^ Among these materials there 
are examples such as A2lr207 (A is a rare-earth element) 
where the relevant effective low-energy degrees of freedom 
are conducting 5d electrons on the Ir^+ ions that live on the 
geometrically frustrated pyrochlore lattice formed by corner- 
sharing tetrahedra.^^"^^ Interestingly, an anomalous Hall ef- 
fect has been observed in the metallic pyrochlore. 

Here, we connect these research directions and consider a 
tight-binding model describing spin-orbit coupled electrons 
on a quasi-two-dimensional slab of pyrochlore including N 




Figure 1. (Color online) The kagome bilayer (a) and trilayer (b). 
The layers are colored differently for clarity and the considered hop- 
ping processes are indicated by arrows in (a). A 2D projection of the 
trilayer is shown in (c) and illustrates the general structure of the py- 
rochlore lattice: the projection has a three-fold periodicity (kagome 
layer m + 3n has the same projection as layer m). 



parallel kagome layers and N — 1 intermediate triangular lay- 
ers (Fig. 1). Our main result is that this model accommodates 
virtually flat bands carrying Chern number C = N that are 
well- separated from all other bands, even for relatively large 
~ 10 (cf. Fig. 2), and as such, provides an intriguing 
platform for yet unexplored phases of matter. 

Setup. Our starting point is the following highly ideal- 
ized model describing Rashba spin-orbit coupled electrons 
on pyrochlore slabs including N kagome layers, /Crn,^ = 
l,...,Ar: 

H ^^Aa^jcT +^ X B.ij)-ao,^c\^Cj^ , (1) 
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Figure 2. (Color online) Bulk dispersion, Es(k), s = 1, 2, . . . , 4A^ — 1, including topological flat bands in the presence of spin-orbit 
coupling, (a) shows a high- symmetry path through the Brillouin zone along which (b)-(i) show the dispersion of (1) for a system with 
N = 1, 2, 3, 4, 8, 12, 50, 100 stacked kagome layers respectively. In each case there is a very flat band (bold orange line) with Chern number 
C = N. The parameters are chosen as ti = — 1, t2 = Ai = 0.3, A2 = 0.2, t± = 1.3. As described in the text, the value of t± is not affecting 
the flat band but can be tuned to move the other bands (thin blue). 



where c]^ creates an electron with spin a at site i and a = 
((Ji, (72, 0-3) is a vector containing the Pauli matrices as its el- 
ements. In the following we consider hopping, tij, between 
(next) nearest neighbors within each kagome plane with am- 
plitude ti (^2) as well as nearest-neighbor hopping, t_\_, when 
involving the triangular layers [cf. Fig. 1(a)]. Moreover, we 
consider only electrons whose spin is aligned in the z direc- 
tion and consider the electric field E^j originating from ions at 
the center of each hexagon in the kagome layer, JCm, to effect 
only the hopping along the (directed) bonds, R^j, within JCm, 
such that = A^jlE^j x R^j| parametrizes the nearest 

and next-nearest spin-orbit hopping strengths Ai and A2. For 
a single layer, = 1, our model (1) reduces to the kagome 
model introduced in Ref. 3. 

In reciprocal space (1) is represented by a (4iV — 1) x 
(4A^ — 1) Hermitian matrix, Hk, which is diagonalized, 
HklV^.(k)) = E,{k)\Mk)), by the states |V^,(k)) = 
J]] • a| (k) where \i) labels the states in the unit cell, s labels 
the bands in order of increasing energies, £^5+1 (k) > £^s(k). 



and i^s = 1,2,..., AN — 1. For an explicit construction of 

we refer to the Supplemental Material. To characterize the 

bands with dispersion £^s(k) we calculate the Chern number, 

^ ~ ^ /bz ^i2(k)<^^^, which is a topologically protected 

integer quantity defined for an isolated band described by the 

wave functions |?/^s(k)), via the Berry curvature, FfjCk) = 

dk^AjCk) — dkjAfCk), which in turn is defined in terms of 

the Berry connection A^(k) = —i{iljs{k.)\dkj\il^sO^))' The 

physical significance of the Chern number is that it counts the 

number of current carrying edge states, and as such gives the 

quantized Hall conductivity of a filled band, an = Cj^^^ 
2 

(and typically cth = Cu^ for an incompressible state at frac- 
tional band filling, u). 

Band structure and wave functions. We begin our discus- 
sion of the band structure of (1) by considering a few instruc- 
tive limiting cases. First, for ti = ±t±^t2 = Ai = A2 = 0, 
there are 2N — 1 perfectly flat bands which can be understood 
from counting the number of linearly independent localized 
modes — for N = 1 these can be taken as localized on each 
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hexagon with amplitudes of alternating sign. (In the three- 
dimensional limit, N ^ oo, this is consistent with the two flat 
3D bands known from bulk pyrochlore^"^.) These perfectly flat 
bands are, however, not carrying a well-defined Chern num- 
ber as dispersive bands touch the flat ones at quadratic touch- 
ing points occurring at k = 0. Next, for = and finite 
spin-orbit coupling. A, we get N decoupled copies of the band 
structure, including C = ±1 bands, discovered in Ref. 3 (in 
addition there are A^ — 1 trivial bands due to the inert triangular 
layers). 

Turning to the full spin-orbit coupled model with tj_ ^ 0, 
we quite generically find a band with Chern number C = N. 
Figure 2 provides examples of the band structure, including 
flat bands with C = (in this example situated ats = N) for 
various N given a suitable choice of hopping parameters. The 
C = N bands have identical (weak) dispersion, Es=nO^), re- 
gardless of the number of kagome layers, N, as will become 
clear later when we consider the properties of the wave func- 
tions. 

To quantify how well the nontrivial bands are separated 
from the other bands we define the flatness ratio, A/W, in 
terms of the bandwidth, W = max [Es(k) — Es{k')], and 



the energy gap, A = min 

k,k'GBZ 



k,k'GBZ 



L(kO,^.+i(kO 



£^s(k)]. A negative flatness ratio does not imply the existence 
of touching points, but it signals that there is no chemical po- 
tential such that only the topological band is partly filled (at 
the noninteracting level). While the topological band itself is 
independent of the value of t± as long as it is nonzero, the 
other bands depend in general on the detailed value of t±, 
leading to a nontrivial -dependence of A/W as shown in 
Fig. 3(a). In Fig. 3(b) we fix = 1.3, ti = -1,^2 = 
0.3, N = 2 and plot A/W as a function of Ai and A2. This 
further illustrates the fact that although reaching optimal val- 
ues is a matter of fine-tuning, there are sizable parameter re- 
gions for which A/W 1. The separation of energy scales 
demonstrated in Figs. 2,3 implies that there is a region of 
parameter space up to relatively large N ~ 10, such that 
A U W, for some interaction strength, U. Note how- 
ever, that one can in principle always make the bandwidth, 
W, arbitrarily small by including longer range hopping^, thus 
the limit on N is not a fundamental one, but rather an esti- 
mate on what might be achieved for reasonable fine-tuning. 
Moreover, adding a chemical potential on the triangular lay- 
ers, H ^ H -\- /jj Xli^K:^ provides additional freedom of 
shifting the other bands while leaving the flat C = N band 
completely unaltered. As in the C = 1 bands studied ear- 
lier, this clearly provides a promising platform for strongly 
correlated states. In particular, the C = > 1 bands could 
lead to new phenomena with no analog in continuum systems 
(Landau levels). 

A prominent feature of the wave functions in the flat 
topological bands is that they are entirely localized to the 
kagome layers. In fact, we find that the wave functions in 
the C = N bands can be completely understood in terms of 
the states |^^^^(k)) in the nearly flat band with Chern num- 
ber C = 1 in the single-layer kagome model. ^ More specif- 
ically, the C = N states can be written as {i/js^nO^)) = 
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Figure 3. (Color online) The flatness ratio A/W (defined in the text) 
as a function of t± for various N (a). Notably, the flatness reaches 
its maximum for relatively large t±. Here we use ti = — 1, t2 = 
0.3, Ai = 0.348, A2 = 0.172. The kink reflects that the upper bands 
become closer to the topological band than the lower bands. In (b) 
we plot A/W for = 2 kagome layers as a function of Ai and A2 
for = 1.3, ti = — l,t2 = 0.3. The flatness ratio is large for a 
sizable region and, in fact, positive throughout the plotted range. 



E™=i«™(k)l<^^='(k)), where \4>g=H^)) is the C = 1 
State localized to JCm- The coefficients am(k) are uniquely 
determined by local constraints similar to those in localized 
modes appearing in frustrated hopping models without spin- 
orbit interactions (cf., e.g., Ref. 34): The relevant amplitudes 
in the kagome layers obey a sum rule leading to a vanish- 
ing total amplitude for hopping to the triangular layers. In 
terms of the wave function amplitudes, this can be stated 



as a. 



where O^. depends on the unit cell and gauge conventions.^^ In 
terms of a^(k), this requirement leads to a nontrivial proba- 
bility distribution, P^(k) = |a^(k)|2 ex e±2m/5(k)^ of being 
in kagome layer JCm, as is illustrated in Fig. 4 and derived in 
the Supplemental Material. 

That the wave functions are completely localized to the 
kagome layers is indeed reminiscent of the structure of multi- 
layer quantum Hall systems, or any system including N de- 



4 




0.8 



0.2 



0.0 





Pi(k) 

-- P2(k) 
P3(k) 

A(k) 














\ 1 
\ ^ t 

'\ ^ : / 




:\\ *. * f/ 

:\\ *, : // 

: \\ * : / // 

1 Vs. J • /y 1 


1 



1 . r 

: - 


A(k) 

-- P2(k) 
P3(k) 














........ 

:\ \ / 
\ \ : / 


\ : / -' 


\\ ■ 1 / \ 1 : 
1 \ \ : / / \ \ [ t ; 

\ / *^ : / t 


r -M -K r KM r 



-M -K 



K M 




Figure 4. (Color online) The probability distribution, Prn(k), of being in different kagome layers, /Cm, in the C = topological bands 
(t± = 1.3, ti = — 1, t2 = 0.3, Ai = 0.3, A2 = 0.2). In (a)-(c) we display the probability distributions for each layer within the C = 2, 3, 4 
bands respectively through a high symmetry path through the Brillouin zone (cf. Fig. 2(a)). That the wave functions are completely localized 
to the kagome layers is manifested by J2m (^) ~ 1, Vk. The wave function is entirely localized to the top or bottom kagome layer at the 
i^-points. In (d) we display the total weight, (Pm), in each kagome layer for various N. Notably, most weight is carried by the outermost 
layers. 



coupled C = 1 bands for that matter. However, the struc- 
ture of the Hilbert spaces is starkly different — while the inter- 
layer dynamics in more conventional multilayer systems can 
in principle be independent, i.e., the layers are coupled only 
via interactions, the particle motion within different kagome 
layers in the C = N bands is necessarily strongly coupled 
even in the absence of interactions between the layers. A man- 
ifestation of this coupling between layers is that translations in 
reciprocal space within the C = bands effectively amount 
to shifting the kagome layer. In Ref. 22, it was argued that 
bands with this type of symmetry enriched translation struc- 
ture can, in principle, harbor entirely new phenomena such as 
wormholes with non-Abelian statistics. 

The nontrivial translation structure is especially transparent 
in the = 2 case highlighted in Fig. 4(a): At the K-point the 
full weight of the wave function is in the lower layer. Pi (k = 
K) = 1, while it is entirely localized to the upper layer at —K, 
P2(k = —K) = 1, i.e., at the i^T'-point. More generally, the 
role of the kagome layers JCm and JCN-\-i-m are interchanged 
by the transformation k ^ — k. 

For C = > 2 the total weight is not evenly distributed 
over all kagome layers (Fig. 4 (b)-(d)). For instance, as a func- 
tion of k, the probability of being in a certain layer reaches 
unity only in the top (/Cat) and bottom (JCi) layers as a direct 
consequence of Prn(k) = |am(k)p ex e^'^^^^^^\ For large 
N this implies that a considerable fraction of the integrated 
weight, (Pm) = ^BZ /bz ^m(k)(i^/c, is carried by the outer- 



most layers, i.e., by the surface as is illustrated in Fig. 4(d). 
This surface localization is also in sharp contrast to more con- 
ventional multilayer systems. 

Discussion. In this work, we have demonstrated that nearly 
flat bands with Chern number C = N exist in a short-range 
tight-binding model on a quasi-two-dimensional slab of the 
pyrochlore lattice. We have also provided a detailed charac- 
terization of the wave functions within these bands and high- 
lighted their structure, which is qualitatively distinct from sys- 
tems hosting multiple C = 1 bands such as multilayer quan- 
tum Hall systems. At partial band filling it is very likely 
that this model harbors intriguing correlated states. Indeed, 
based on an elegant, albeit rather implicit, Wannier function 
approach^^, it was recently argued^^ that lattice dislocations 
can lead to qualitatively new phenomena that are only real- 
izable in C > 1 bands. Here, we have provided a promising 
platform for realizing such phenomena and a testbed for future 
numerical studies of flat bands with variable Chern number. 

An intriguing open question is whether our model, mutatis 
mutandis, can be realized experimentally. In this context it 
is indeed promising that there is a rapid development in en- 
gineering quasi-two-dimensional (multi-layer) oxides^^, and 
in particular, that there exists a number of pyrochlore-based 
transition-metal oxide materials, especially iridates, where the 
conduction electrons experience strong spin-orbit effects. 

It has recently been shown that our model indeed harbors a 
large number of new fractional Chern insulator phases. 
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SUPPLEMENTARY MATERIAL FOR "FLAT BANDS WITH 
HIGHER CHERN NUMBER IN PYROCHLORE SLABS ' 

In this supplementary material, we provide details on how 
to construct the Hamiltonian in reciprocal space as well as a 
definition of the unit cell. We also derive an explicit expres- 
sion for the A^-layer C = N states in terms of the single-layer 
C = 1 states and explain the exponential localization of the 
wave functions to the surface layers. Finally, we present Berry 
curvature data for various C = N bands. 



Unit cell, gauge conventions and construction of the 
Hamiltonian 

To derive the reciprocal space Hamiltonian, we need to 
fix the gauge and choose an appropriate unit cell. Here we 
chose a unit cell as depicted for A/" = 3 in Fig. 5 which has 
the benefit that the construction of generalizes in a straight 
forward way to any number of kagome layers, N. Further, we 
use the gauge freedom such that hopping inside one unit cell 
does not give rise to a phase factor, while hopping a distance 
d = nai + 1713.2 (to another unit cell) gives rise to a phase 
factor e*^'^. With these conventions, and /c^ = k • a^, z = 
1, 2, /c3 = /ci — /c2, the Hamiltonian for one kagome layer in 
reciprocal space reads 



njN=l 
Atk 



+ zAi 



t2 e- 
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which is related to the form given in Ref. 3 by a gauge trans- 
formation (and a change of basis). 

Our choice of unit cell makes the generalization to mul- 
tiple coupled kagome layers especially simple as it allows 
us to reuse H^^^ for each of the kagome layers. More- 
over, the phases acquired when hopping to the triangular sites, 
z = 4, 8, 12, AN — 4, are the same for each triangular layer 
using these conventions. 

For the two layer system one finds 
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Figure 5. Generalizable unit cell, shown for three layers. The basis 
vectors are ai = 2ax, a.2 = ax — V^ay, where a is the lattice 
constant. 



which can be more compactly by introducing 

/ 1 \ 
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= ti_ 
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-iks 



(4) 



V 1 / 



with which we can write (3) as: 

/Ar=i 



njN=2 
Atk 









iN-- 



(5) 



This generalizes directly to an arbitrary number of layers, N, 
yielding the following (4A^ — 1) x (4A^ — 1) matrix: 
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The eigenstates in band s are 



/ af(k) \ 
a^(k) 

V«liv-i(k)y 



(6) 



(7) 



which can also be written as |^s(k)) = a|(k)|i), as in the 
main text. 
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Wavefunction structure and surface localization 

With the conventions above, the constraints for the C = N 
bands, 

^ e<a^^=^(k) = Vk,m, (8) 

iG{/C^U/C^+i} 

where Km U /Cm+i = {3m — 2, 3m — 1, 3m, 3m + 2, 3m + 
3, 3m + 4} are fulfilled for 



oi = o 



for j = 3m + 2 (9) 
for j = 3m + 3 (10) 
otherwise. (11) 



We can now explicitly express the state |V^s=Ar(k)) = 
Sm=i <^m(k)|0^^^(k)) in terms of the single layer eigen- 
state |^^=^(k)) := (ai(k), a2(k), a3(k))^. Together with 
(8) this leads to the relation 



Q^m+i(k) 

a^(k) " e-^^2ai(k) + e^^3a2(k) +a3(k) 



(ai(k)+a2(k)+a3(k)) 



r(k). 



(12) 



We emphasize that the ratio r(k) = """(i^^^ does not depend 
on TO. From this we directly get 

a™(k) = ao(k)r(k)'" , (13) 

where ao(k) is set by normalization. One finds 

AT AT 



1 = J2 = E l«o(k)r(k) 



m\2 



N-1 



|ao 



(k)|2 E (IKk)n™+^ (14) 



l«o(k)P = 



|r(k) 



|r(k)|2 (|r(k)|2)Af - 1 



(15) 



except for |r(k)p = 1 where = l/A/". Now we can write 
the eigenstates in the C = TV band as 



AT 



\i^s=NO^)) = E ao(k)r(k)™|(/.^=i(k)) . (16) 



i.e., as an entangled sum of the wave functions in the single- 
layer C = 1 bands. 



To highlight the surface localization we define the k- 
dependent localization length, (5(k) = 1/| log(|r(k)|)|, which 
leads to 



P^(k) = |a^(k)|^ = |ao(k)|^e 



2 ±2m/5(k) 



(17) 



as stated in the main text. 

Berry curvature 



In Fig. 6, we display the Berry curvature for various 
N given the conventions specified above. Compared to the 
C = N = 1 band (Fig. 6(a)) the maximum of the Berry cur- 
vature gradually shifts position for increasing TV (Fig. 6(b)- 
(f)) . Meanwhile, it remains smooth, positive (semi-)definite 
and free from singularities. 



8 




Figure 6. (Color online) Berry curvature in the C = bands with ti — — 1, t2 = 0.3, Ai = 0.3, A2 = 0.2 (independent of ti_). In (a)-(f) we 
display the results for = 1, 2, 3, 6 respectively. 



